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Mass Estimation without using MET in early LHC data 
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Many techniques exist to reconstruct New Physics masses from LHC data, though these tend 
to either require high luminosity 0(100 fb~^), or an accurate measurement of missing transverse 
energy (MET) which may not be available in the early running of the LHC. Since in popular models 
such as SUSY a fairly sharp, triangular dilepton invariant mass spectrum can emerge already at low 
luminosity 0(1 fb~^), a Decay Kinematics (DK) technique can be used on events near the dilepton 
mass endpoint to estimate squark, slepton, and neutralino masses without relying on MET. With 
the first 2 fb~^ of 7 TeV LHC data SPSla masses can thus be found to 20% or better accuracy, at 
least several times better than what has been taken to be achievable. 
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We in the particle physics community are naturally ea- 
ger to glimpse signs of New Physics in early data from 
the running of the LHC, most likely to be first seen 
in anomalous values of inclusive measurements of lep- 
ton and hadronic jet activity, accompanied by missing 
energy from an escaping Dark Matter particle. How- 
ever, to glean more quantitative information such as 
new particle masses will require fitting data to assumed 
decay topologies, and a plethora of such mass recon- 
struction techniques have accumulated over the years [li]. 
Typically, these require either large datasets, with inte- 
grated luminosities f £ ~ 0(100 fb~ ), as in endpoint 
formula techniques [2|, |3|, or precision measurement of 
missing transverse energy/momentum (MET), as in the 
mT2[4-6] and polynomial methods [7|-|9|. However, accu- 
rate MET measurements are not expected to be available 
in early LHC datailQj. 

The Decay Kinematics (DK) technique, characterized 
by the full reconstruction of events that lie near the end- 
point of an invariant mass distribution for the decay 
products and where the kinematics are exactly known, 
has made its debut recently and proven useful as a mass 
estimator in the commonly considered scenario of neu- 
tralino pair production[l llll2| . assuming accurately mea- 
sured MET at high luminosity. The purpose of this Let- 
ter is to demonstrate that DK may in fact be the tech- 
nique of choice in early analysis of certain decay chan- 
nels, in particular the well-studied squark-initiated cas- 
cade that ends in the lightest supersymmetric particle 
(LSP) and Dark Matter candidate, the lightest neutralino 
X?: 
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We will show that the DK technique can outperform cur- 
rent methods when one can rely neither on MET nor large 
datasets. 

Consider a collection of N events, each having at least 
one pair of isolated opposite-sign same-flavor (OSSF) lep- 
tons plus a high energy jet {i^£^ j), which we will assume 



to have arisen from the aforementioned squark cascade 
([1]). For each event, the following mass-shell constraints 
hold in the narrow-width approximation: 
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where P^^, P^f, and ^n are the four- momenta in the n"^ 
event (n = l,...,iV) of the jet, leptons and the LSP, 
and where rrig, TO2,wi and rus are abbreviations for the 
relevant squark, neutralino, and slepton masses. Note the 
ambiguity in which lepton is assigned to which sparticle 
decay. We will return to this below. 

For events that have exactly maximum dilepton invari- 
ant mass, i.e. if 

{P,t + Pn? - Mu ^ {ml - ml){ml - ml)/ml (6) 

then as shown in [12| the longitudinal component of the 
LSP's three-momentum is subject to a coplanarity con- 
straint: 



/n^ — ^ -y nx i y p^ny , 
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In the spirit of [13], for each event the constraints (I1])-(IS]) 
can be rearranged into three equations linear in the four 
components of _^„ : 



{P^+Pn+Pn? + 2 A • [Pi+P^^ + P-) = A'-l, (9) 

{P++p-f + 2 A • (F+ + P-) = A^i, (10) 

2pr.-P- ^ Al,. (11) 



with At =ml-m\, All 



ml, and A^^ 



These, when combined with ([7]), can be used to 



solve for pn given a mass hypothesis {mq,m2,ms,mi). 
Thus, for a collection of perfect events arising from ([I}, 
the quantity 



(12) 



would be exactly zero if the mass hypothesis were correct. 
Note that in this procedure the ambiguity in identifying 
which lepton comes from the slepton decay can be re- 
moved by explicitly reconstructing the X2 ^^st frame via 
the technology of (I4I and seeing which lepton is paral- 
lel/antiparallel to the LSP. We will use this below. 

Of course, in a real data sample no event will lie ex- 
actly at the kinematic endpoint, and one must settle for 
a collection of events within some window of size e near 
the dilepton maximum M^. Experimental effects in the 
measured momenta will introduce further smearing on 
the solution, but it is reasonable to believe that the cor- 
rect neighborhood in (niq, 1712, rris, mi)-space should still 
be about that point which minimizes ^. 

The search over this mass space can be greatly simpli- 
fied by assuming ((6]) to hold exactly and using informa- 
tion from estimates of the upper endpoints of the qi£ and 
qi invariant mass distributions, Mqu and Mqi, to indi- 
cate the kinematically allowed regions, as theseprovide 
analytical constraints on the unknown masses [2!, l3|. Even 
a very crude estimate of these endpoints (call them Mqu 
and Mqi) will be sufficient input for the DK technique, 
as we will see below. 

We have tested this method in a realistic Monte Carlo 
simulation of the SUSY benchmark point SPSla[lJ| at 
low integrated luminosity (2 fb~ ) and center-of-mass en- 
ergy 7 TeV, as expected from the first first two years of 
LHC running. For both signal and background produc- 
tion we use PYTHIA 6.413 [15J interfaced to the fast sim- 
ulation of a generic LHC detector, AcerDET-1.0 16j. For 
further details of the simulation, see 1^ . The cuts used 
to isolate our signal are: each event must have at least 
two hard jets with pT > 150, 100 GeV (expected from 
pair production of squarks or gluinos), plus two isolated 



OSSF leptons (e+e" or fi+n') with pr > 10 GeV. We 
apply pt dependent lepton efficiencies as in[l2|, based on 
full simulation results published in[17|. The jet resolution 
at 7 TeV is not yet well understood, and one possible 
parametrization of this is to introduce additional smear- 
ing of jet momenta by hand; yet the uncertainty in jet 
direction is likely to be much smaller than of jet energy, 
so it is perhaps more realistic to scale all components of 
jet 4-momenta by the same factor — we thus scale all jet 
energies up by a generous 5% (scaling down gives virtu- 
ally the same influence on final results). Standard Model 
backgrounds for the jet-|-dilepton endstate considered are 
overwhelmingly dominated by tt which we have gener- 
ated for this particular machine energy and luminosity 
(PYTHIA (7LO = 85 pb), giving 1.7 • 10^ background 
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FIG. 1: Events from the shaded bins of the OSSF distribution 
(74 GeV < Mu < 78 GeV) are used for DK analysis. 



events; lesser backgrounds include ZZ* and Z + jets 
which are several times smaller and can be essentially 
eliminated by selecting events away from the Z-pole. All 
of this is against only 1.2-10^ SUSY events (all processes), 
but we will nevertheless see below that DK has excellent 
background rejection. 

As can be verified from Figll] even at this low luminos- 
ity and set of cuts the OSSF dilepton distribution typ- 
ically contains hundreds of events (here about 600, half 
of which are from SUSY), and dozens of events within a 
moderately small window (here chosen to be e = 4 GeV) 
near the kinematic endpoint at Mu « 77 GeV. Note that 
although backgrounds outnumber signal in general, this 
is not so for events near the dilepton endpoint: in this 
sampling area shown with 74 GeV < Mu < 78 GeV, 24 
of the events are from SUSY while only 4 are background 
(S/B = 6), hence choosing events near this endpoint gives 
us quite some leverage against backgrounds. In fact, ow- 
ing to the fact that background events tend to fail out of 
the DK algorithm (see below), we find that backgrounds 
only become a serious issue when S/B ~ 1 and S < 10. 
We must mention that one has some freedom in choosing 
e, balancing the tension between accurate results with 
small systematic errors (small e) and low statistical error 
(larger e). As a rough criterion in general, one should 
choose e so as to give C'(10)-C'(100) events; our experi- 
ence is that this suffices to give a decent mass estimate. 

To determine what set of SUSY masses is most con- 
sistent with this sample of events, we perform the fol- 
lowing scan over mass space: rus is chosen in a liberal 
range, 100 GeV < nis < 500 GeV (lower bound from 
LEP constraints) as well as Mqu and Mqi in ranges 
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FIG. 2: Mqi and Mqti mass distributions, flavor-subtracted as 
described in the text and in Q] • A rough estimate of the end- 
points suffices for our analysis: 330 GeV < Mqu < 430 GeV, 
300 GeV < Hal < 400 GeV. 



estimated from FigH 300 GeV < M^i < 400 GeV, 
330 GeV < Mqu < 430 GeV. Note that, as is cus- 
tomary in handling jet ambiguity, what we have actu- 
ally plotted are the minimum values of Mqu and Mqg 
computed from each of the two highest energy jets; also, 
before this minimization Alqi is taken to be the larger 
of the two possible q£ masses. To reduce background 
and combinatorial effects further, these distributions are 
also computed for OSOF lepton events, labeled Mqw and 
Mqii, and are subtracted from the Mqu and Mqi distri- 
butions, respectively. The total procedure thus improves 
endpoint precision at the cost of distorting distribution 
shape somewhat, but this is irrelevant to our present pur- 
pose. 

For each choice of the endpoints M«, Mqu and Mqi, 
the other unknown masses are then fixed from analytic 



expressions of the endpoints in terms of these masses and 
TOs as 



mi 



m,Mu\/Mqff/M^/, 



m2 = m^Y l-f M^^/(m2 - mf), 



rriq = m2^Jl + M u/{m 
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(13) 

(14) 
(15) 



and for each such point in mass space the linear system 
of equations ([7l)-(fTT|) is solved event- wise for /P„, con- 
structing the sum (fT2l) defining ^. Both possible lepton 
assignments are tried. In this procedure each event con- 



tributes the minimum value of 



-mi I resulting from 



matching the OSSF lepton pair to each of the two highest 
energy jets, accepting only solutions where the assumed 
lepton identity does indeed lead to reconstruction of the 
correct leptonic directions relative to the LSP momen- 
tum, i.e. the cosine of the angle between the assumed 
near(far) lepton momentum and the LSP must be >(<)0 
in the rest frame of the decaying X2- I^: fo^ ^^Y event 



the quantity 



mi I exceeds a tolerance T (we 



take r = 30 GeV), or if /P„ is tachyonic, or if the re- 
construction of the X2 velocity fails for both possible jet 
pairings, the event is assigned a fixed contribution of T 
to ^ {i.e. a fit which is T = 30 GeV in error is deemed as 
bad as finding no fit at all). As with the e parameter, T 
is adjustable and could be optimized somewhat, but any 
reasonable value will do. The choice of {m,s,Mqu,Mqe) 
giving the minimum ^ is then assumed to yield the 
correct mass hypothesis. For the particular data set 
shown, we obtain {ms,Mqee,JIqi) = (133,408,300) GeV 
as the ^-minimizing solution, and using (fT3l) - (fT5l) yields 
{mi,ms,m2,mq) = (90.3,133.0,170.0,510.7) GeV, quite 
close to the nominal values of (96, 143, 177, 537) GeV. 

To get some idea on the expected accuracy of this 
method owing to random fluctuations in signal and 
background, we look at the best-fit solutions in 50 
independently-generated data sets, obtaining the follow- 
ing results: 

mi = 89.9 ±21.8 GeV (16) 

m, = 131.0 ±22.2 GeV 

ma = 169.2 ±21.6 GeV 

mq = 512.4 ±31.4 GeV 

Quoted uncertainties are actually square-roots of sample 
variances. Though these results are slightly lower than 
the nominal values, one generally gets accurate (within 
^ 15 — 20 %) estimates of the masses, and the spread 
in the solutions is consistent with the rough width of the 
peak in, for example, the LSP mass distribution in Fig. [S] 
This width scales with the jet energy-bias we put in by 
hand, so once low-luminosity jet reconstruction efficiency 
is better understood the accuracy of the method would 
be closer to '^ 10 — 15% on average. 
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FIG. 3: Best fits to the LSP mass for each of 50 
independently-generated 2 fb~^ data sets (SUSY signal + SM 
background). 



only that one can utilize the high statistics in the dilep- 
ton endpoint, but that one can resolve lepton ambiguity 
through explicit reconstruction that does not rely on a 
measurement of MET. In fact, the same method may be 
used on a decay chain where the LSP decays through R- 
parity violating operators, as there is no reliance on pair- 
production of new sparticles. The need for an effective 
mass estimator in the early stages of the LHC is attrac- 
tive not only to help theorists hone in on likely models for 
New Physics signals, but also in the design and optimiza- 
tion of experiments, especially with respects to the LSP, 
whose mass measurement at the LHC will directly affect 
designs of Dark Matter detection facilities elsewhere. 

Thanks to S. Kraml for some useful suggestions, and 
to A. Raklev for helping with background data and giv- 
ing much invaluable critique of the manuscript. MJW is 
supported by the Australian Research Council. 



Let us again stress that, prior to this Letter, the only 
known method of reconstructing particle masses in early- 
LHC data was the standard endpoint analysis, and this 
will give considerably worse results: in addition to the 
three endpoints M«, Mqu, and Mgi whose precise de- 
termination is not possible with such low statistics, one 
needs a fourth endpoint, usually taken as the lower end- 
point of the Mqii distribution, to solve for a point in 
(rnq,m2,mg,mi)-spa,ce. However, this lower endpoint, 
which is supposed to be at Mgu = 202 GeV, is not dis- 
cernible in Fig. [2] Studies with more sophisticated detec- 
tor simulation and endpoint fitting algorithms point to 
uncertainties in the masses that can reach 50% or more 
from early LHC data[10J at comparable luminosities. 

This Letter demonstrates that, in the early stages of 
the LHC, one may obtain far more accurate estimates 
of New Physics masses than previously thought if one 
uses a DK technique, at least for i~^i~j events assumed 
to arise from the squark cascade chain considered here. 
Generalization to the same final state arising from other 
models (such as extra dimensions or little higgs mod- 
els), or to different decay topologies is straightforward. 
What makes DK so powerful in the present case is not 
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